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ABSTRACT 

It is widely suspected that AGN activity ultimately sweeps galaxies clear of their gas. We work 
out the observable properties required to achieve this. Large-scale AGN-driven outflows should have 
kinetic luminosities ~ r/LEdd/2 ~ 0.05LEdd and momentum rates ~ 20LEdd/c, where I/Edd is the 
Eddington luminosity of the central black hole and rj ~ 0.1 its radiative accretion efficiency. This 
creates an expanding two-phase medium in which molecular species coexist with hot gas, which can 
persist after the central AGN has switched off. This picture predicts outflow velocities ~ 1000 — 
1500 kms -1 and mass outflow rates up to 4000 Moyr -1 on kpc scales, fixed mainly by the host 
galaxy velocity dispersion (or equivalently black hole mass). All these features agree with those of 
outflows observed in galaxies such as Mrk231. This strongly suggests that AGN activity is what 
sweeps galaxies clear of their gas on a dynamical timescale and makes them red and dead. We suggest 
future observational tests of this picture. 

Subject headings: galaxies: evolution — quasars: general — black hole physics — accretion 



1. INTRODUCTION 

Recently, thr ee grou ps (iFeruglio et al. 120101 : 

IRunke fc Veilieu^ 12011 ISturm et all 1201 lb have 
used molecular spectral line observations to re- 
veal fast (v out ~ 1000 km s" 1 ) kpc-scale, massive 
(M out ~ 1000 Moyr^ 1 ) outflows in the nearby quasar 
Mrk231. Other galaxies show indications of simil ar 
phenomena (e.g. iRiffel fc Storc hi-Bcrgmannl l2011bl lal: 
ISturm et al.l l201lh . These appear to show how quasar 
feedback can transform young, star-forming galaxies 
into red and dead spheroids. All three groups reach 
this conclusion for Mrk231 essentially by noting that 
the mass outflow rate M out and the kinetic energy 
rate E out = M out vl ut /2 of the outflow are too large 
to be driven by star formation, but comparable with 
those predicted in numerical simulations of AGN 
feedback. The kinetic energy rate is a few percent 
of the likely Eddington luminosity iEdd = 4:TtGMc/k 
of the central black hole, of mass M (where k is the 
electron-scattering opacity). The outflowing material 
must have a multi-phase structure, because v ou t greatly 
exceeds the velocity corresponding to the molecular 
dissociation temperature (i>di ss ^ 10 km s , see Section 
5 below). 

In a recent paper ([King et al.l 120111) we showed that 
large-scale flows of this type (technically, an energy- 
driven flow, see Section 3) can indeed drive much of the 
interstellar gas out of a galaxy bulge on a dynamical 
timescale ~ 10 8 yr, leaving it red and dead, provided that 
the central supermassive black hole accretes for about 
twice the Salpeter time after reaching the value set by 
the M - a relation. In iPower et all fl2011l) we showed 
that the remaining bulge mass is close to the value set 
by the observed bl ack-hole - bulge-mass relation (e.g. 
IHaring & Rixll2004[ ). However we did not investigate the 
observable features of this process, including in particu- 
lar the way that the interstellar gas is swept up. 
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We return to this problem here, as it offers a clear 
observational test of the idea that AGN outflows are re- 
sponsible for making galaxies red and dead. To keep our 
treatment as general as possible (specifically, indepen- 
dent of the details of numerical simulations) we adopt 
a simple analytic approach. We find that this process 
predicts outflow velocities ~ 1000 — 1500 kms" 1 , and 
mass outflow rates up to ~ 4000 M Q yr -1 , several hun- 
dred times the Eddington value, in good agreement with 
observations. In addition, we find that the observable 
momentum outflow rate is ~ 20 times greater than L/c 
of the driving AGN, also in agreement with observations. 
We conclude that AGN outflows are good candidates for 
the agency sweeping galaxies clear of gas. 

2. WINDS 

To drive outflows with E out approaching L^dd at large 
radius, the active nucleus of a galaxy must somehow 
communicate this luminosity from its immediate vicin- 
ity. Direct transport by radiation is problematic, not 
least because galaxies are generally optically thin (How- 
ever, dust opacity may be large e nough to absorb thi s 
radiation and provide feedback; see iMurrav et al.l feOOS') . 
Jets are sometimes invoked, but are relatively inefficient 
because they tend to drill holes in the interstellar medium 
rather than driving it bodily away. Accordingly the most 
likely connection is via high-velocity wide-angle winds 
expelled from the vicinity of the nucleus by radiation 
pressure (e.g. iPounds et al.l 2003a. b). Recent observa- 
tions suggest that suc h winds are very common in AGN 
(jTombesi et al.H2010allH ). In this paper, we use the term 
'wind' to refer to the mildly relativistic (v ~ 0.1c) ejec- 
tion of accretion disc gas from the immediate vicinity 
of the SMBH resulting from Eddington accretion, and 
'outflow' (see Section 4) for the large-scale nonrelativis- 
tic flows caused by the interaction between the wind and 
the galaxy's ambient gas. 

The winds have simple properties. With mass rate 
M w ~ MEdd, where Msdd = ^EddA/c 2 is the Ed- 
dington accretion rate and rj is the accretion efficiency, 
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Fig. 1. — Diagram of momentum-driven (top) and energy-driven (bottom) outflows. In both cases a fast wind (velocity ~ r\c ~ 0.1c) 
impacts the interstellar gas of the host galaxy, producing an inner reverse shock slowing the wind, and an outer forward shock accelerating 
the swept— up gas. In the momentum-driven case, the shocks are very narrow and rapidly cool to become effectively isothermal. Only the 
ram pressure is communicated to the outflow, leading to very low kinetic energy ~ (c/c)Le ( j ( j. In an energy— driven outflow, the shocked 
regions are much wider and do not cool. They expand adiabatically, communicating most of the kinetic energy of the wind to the outflow 
(in simple cases approximately l/3rd is retained by the shocked wind). The outflow radial momentum flux is therefore greate r than that 
of the wind. Momentum-driven flows occur when shocks happen within ~ 1 kpc of the AGN, and establish the M — cr relation (King 2003, 
2005). Once the supermassive black hole mass attains the critical M — a value, the shocks move further from the AGN and the outflow 
becomes energy— driven. This produces the observed large— scale flows, which probably sweep the galaxy clear of gas. 



a wide-angle wind h as scattering optical depth ~ 1 
(jKing fc Pounds! 12001 . assuming that the covering fac- 
tor of absorbing gas is close to unity (see the discussion 
below eqn[5]). So each driving photon on average scatters 
about once before escaping to infinity and gives up all of 
its momentum to the wind, so that the wind mass flow 
rate M w and velocity v obey 
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Defining m = M w /MEdd ~ 1 as the Eddington factor of 
the wind, we immediately find 



i ~0.1 

c m 



(2) 



(cf lKingll2010l ). The likely ionization e quilibrium o f the 
wind is such that it produces X-rays (jKingl [20101 ) . In 
line with these expectations, blueshifted X-ray iron ab- 
sorption lines corresponding to velocities ~ 0.1c are seen 
in a si gn ificant fraction of loca l AGN (e.g. iPounds et al.l 
I2003aia [Tombesi et al.|[2010allH ) . justifying our assump- 
tion of a covering factor close to unity. In all cases the 
inferred wind mass flow rates agree with eqn. ([T]). So 
these black hole winds have momentum and energy rates 
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where we have used eqns |TJ [5J) in eqn ([J} . 

3. SHOCKS 

The expression (Q| for the energy rate of a black hole 
wind is obviously promising for driving the observed 
large-scale outflows. Although the interstellar medium 
is clumpy, the outflow bubble inflated by the wind eas- 
ily sweeps past the clumps , affecting the diffuse gas (e.g. 
iMac Low fc McCravlll988() . Furthermore, the clouds are 
shocked by the p assing outflow and eva porate inside the 
hot wind bubble (jCowie fc McKedll977t) . so most of their 
material also joins the outflow. A detailed treatment of 
the interaction between the clumpy ISM and the wind- 
driven outflow is beyond the scope of this paper, although 
we address some of the implications in Section 6. In the 
present analysis we assume that most of the sightlines 
from the SMBH are covered with diffuse medium, irre- 
spective of whether they are also obscured by clumps. 

The question now is how efficiently the wind energy 
is transmitted to the outflow. This depends crucially 
on how the wind interacts with the diffuse interstellar 
medium of the host galaxy. Since the wind is hyper- 
sonic, it must decelerate violently in a reverse shock, and 
simultaneously drive a forward shock into the host ISM. 
There are two possible outcomes, which are realised un- 
der different conditions in galaxies. 

The first outcome (momentum-driven flow) occurs if 
the shocked wind gas can cool on a timescale short com- 
pared with the motion of the shock pattern. In this 
case the shocked wind gas is compressed to high density 
and radiates away almost all of the wind kinetic energy 
(i.e. i?out << E w — (r/ / 2m) LEdd) ■ This shocked wind 
has gas pressure equal to the pre-shock ram pressure 
P w ~ LEdd^/c oc M, and this pushes into the host ISM. 

The second case (energy-driven flow) occurs if the 
shocked wind gas is not efficiently cooled, and instead 
expands as a hot bubble. Then the flow is essentially 
adiabatic, and has the wind energy rate, i.e. i? ut — 
= (?7/2m)L Edd ~ 0.05L Edd (from eqnHJ). The hot 
bubble's thermal expansion makes the driving into the 
host ISM more vigorous than in the momentum-driven 
case. Observed galaxy-wide molecular outflows must be 
energy-driven, as demonstrated directly by their kinetic 
energy content (cf eqn|4]). 

Which of these two very different cases occurs at a 
given point depends on the cooling of the shocked gas. 
It is easy to show that the usual atomic cooling pro- 
cesses (free-free and free-bound radiation) are negligible 
in all cases. The dominant process tending to cool the 
shocked black hole wind is the inverse Compton effect 
(jCiotti fc Ostriked [19971) . The quasar radiation field is 
much cooler than the wind shock temperature (typically 
~ 10 7 K and ~ 10 11 K respectively), and so cools the 
shocked wind provided that it is not too diluted by dis- 
tance. Equations 8 and 9 of iKind (|2003l) show that this 
holds if and only if the shock is at distances R < 1 kpc 
from the active nucleus, since the Compton cooling time 
goes as R 2 and the flow time typically as R. So we ex- 
pect momentum-driven flow close to the nucleus, and 
energy-driven flo w if gas can be d riven far away from it. 

We note that iSilk fc Nusserl (|2010T) claim that an 
energy-driven flow never occurs. However they seem to 
have considered the cooling of the ambient gas, rather 



than the shocked wind and swept-up ISM: the cooling 
function they use extends only to temperatures ~ 10 7 K, 
rather than the w ind shock temperature ~ 10 11 K - see 
IKing et all (|2011[ ) for details. 

4. OUTFLOWS 

For an isothermal ISM density distribution with ve- 
locity dispersion a and gas fraction f c (the ratio of gas 
density to background potential density) one can solve 
analytically the equation of motio n for the sh o ck pa ttern 
for both momentu m-driven flow (|Kin g 2003, 2005]) and 
energy-driven flow (|King|l2005l; IKing et al.ll2011l ). 

In the momentum-driven case there are two distinct 
flow patterns, depending on the black hole mass M. For 
M < M,, where 
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(with f c = 0.16 (its cosmological value) and CT200 — 
cr/(200 kms -1 )) the wind momentum is too weak to 
drive away the swept-up ISM, and the flow stalls at some 
point. For M > M a the wind momentum drives the 
swept-up matter far from the nucleus. It is intuitively 
reasonable to assume that the black hole cannot easily 
grow its mass significantly beyond the point where it ex- 
pels the local interstellar gas in this way, i.e. beyond M a . 
Equation ([5]) is very close to the observed M — a relation, 
despite having no free parameter. Detailed calculations 
(jPower et al.ll201lD show that the SMBH is likely to grow 
for 1-2 additional Salpeter times after it reaches M a , in- 
creasing its final mass by a factor of a few. This process 
is even more pronounced at higher redshift, as then it 
takes longer for the outflow to clear the galaxy, so the 
SMBH must be active for longer. 

We conclude that outflows drive gas far from the nu- 
cleus, and thus become energy-driven, once M > M a . 
This is evidently the case needed to explain the molecu- 
lar outflows seen in Mrk231 and other galaxies. 

5. LARGE-SCALE FLOWS 

In an energy-driven flow the adiabatic expansion of the 
shocked wind push es the swept-up interstellar medium 
in a 'snowplow'. IKing et all (|201lD derive the analytic 
solution for the expansion of the shocked wind in a galaxy 
bulge with an isothermal mass distribution. With AGN 
luminosity Z^Edd, all such solutions tend to an attractor 
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925a 2 2 / 3 (lf c /f g ) 1/3 km s" 1 
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until the central AGN luminosity decreases significantly 
at some radius R = Ro, when the expansion speed decays 
as 

where x = R/Ro > 1- In eq. ((6]), fg is the gas fraction 
relative to all matter. This may be lower than the value 
f c prevailing when the earlier momentum-driven outflow 
establishes the M — a relation ([5]), as gas may be depleted 
through star formation for example. 

The solutions [7]) describe the motion of the contact 
discontinuity where the shocked wind encounters swept- 
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TABLE 1 
Observed outflow parameters 



Ubject 


i\,r /Ayr 


<t/ km s * 


Lbol/ erg s 1 (/) 


-Mout/MQ yr 


v out /km s 1 


Mrk23l( a > 


4.7- lO 7 ' 6 ' 


120^) 


45.69< c ) (0.80) 


420 


1100 


Mrk23l( d ) 


4.7- 10 7 


120 


45.69 (0.80) 


700 


750 


Mrk23l( e > 


4.7- 10 7 


120 


46.04^' (1.8) 


1200 


1200 


IRAS 08572+3915^ 


~ 4.5 ■ 10 7 * 


120*" 


45.66 (1*) 


970 


1260 


IRAS 13120-5453< e > 


5.3 ■ 10 6 * 


70*** 


44.83 (1*) 


130 


860 


IRAS 17208-0014( e )** 






45.11 « 1) 


90 


370 


Mrkll57(f) 


8.3 ■ 10 6 


100 


42.57 (3.4 ■ 10" 3 ) 


6 


350 


2QZJ002830.4-281706( h ) 


5.1 • 10 9 W 


385*** 


46.58 (5.8 ■ 10" 2 ) 


2000 


2000 



Outflows observed in molecular gas (Mrk231, IRAS 08572+3915, IRAS 13120-5453 and IRAS 17208-0014), and warm ionised gas 
(Mrkll57, 2QZJ002830. 4-281706). * - the A GN is assumed to be radiating at its Eddington limit; ** - the galaxy is starburst- 
dominatcd (Riffcl & Storchi-Bcrgmann 2011b), so we expect a low Eddington factor and henc e make no estimates; * ** - the SMBH 
is assumed to lie on the M — a relation; note that this may be questionable in some cases (cf McConncll ct al. 2011). 

References: a - IRupke fc Veilleuxl 120TH ); b - ITacconi eFaTI 120021) ; c - ILonsdafe et~aH 12003) ; d - IFeruglio eFaLI 120101) ; e - ISturm et at] 
120111) ; f - 1 Veilleux et~aT1 pM) ; 9 - IRiffel fc Storchi-Bergrnamil l !2011bl) ; h - ICano-Diaz et al.l 120111 ); 1 - IShemmer eFaLl 1200? ). 

TABLE 2 

Observationally derived versus theoretically predicted outflow parameters 



Object 


Bout 

0.05L boI 


M out v out c 
Lbol 


fh = 


Mout 
M acc 


/i^pred. 


A/prcd./M yr" 1 


Upred./ km s 1 


Mrk231 


0.66 


18 


490 = 


= 22 2 


840 


880 


810 


Mrk231 


0.51 


20 


820 = 


= 29 2 


840 


880 


810 


Mrk231 


1.0 


25 


1400 


= 37 2 


1110 


1150 


1060 


IRAS 08572+3915 


2.1 


50 


1200 


= 35 2 


910 


950 


875 


IRAS 13120-5453 


0.88 


31 


1080 


= 33 2 


1870 


220 


610 


IRAS 17208-0014 


0.06 


4.9 


396 = 


= 20 2 








Mrkll57 


1.3 


110 


9270 


= 96 2 


170 


85 


115 


2QZJ002830.4-281706 


1.3 


20 


307 = 


17.5 2 


74 


8200 


740 



The first three columns give quantities derived from observations of large-scale outflows summarized in Table [T1 T h e la st three 
columns give the mass-loading parameter, mass sweep— out rate and terminal velocity predicted by our equations ||14|I . l|15|l and (Jsj 
respectively. We assume the simplest case of an isotropic outflow. Collimation would reduce the predicted mass outflow rate and 
increase the predicted outflow velocity. With one outlier (see below), the outflow kinetic energy is always very close to 5% of L^oi 
(1st column) as predicted by eq. l(4}, and the mo men tum loading (2nd column) is always very similar to the square root of the mass 
loading (rhs of 3rd column), as predicted by eq. 1181 1. It is striking that the relation holds for local quasars (Mrk231), high-redshift 
quasars (2QZJ002830. 4-281706) and low luminosity galaxies (Mrkll57). The last two columns can be directly compared with the 
last two columns of Table[T] the discrepancies arise due to strong outflow collimation. The only significant outlier, IRAS 17208-0014, 
is known to be a starburst-dominated galaxy, so we would not expect the outflow to be dominated by the AGN contribution. 



up interstellar gas (see Figure [lj . The observed molecu- 
lar lines are likely to come from the shocked interstellar 
gas ahead of this discontinuity - its temperature is much 
lower (~ 10 7 K) than that of the shocked wind, as we 
shall see. The outer shock must run ahead of the con- 
tact discontinuity into the ambient interstellar medium 
in such a way that the velocity jump across it is a factor 
(7 + l)/(7 — 1) (where 7 is the specific heat ratio). This 
fixes its velocity as 



Vout 



7 + 1 



R ~ 1230cr 



2/3 
200 



1/3 



km s 



(8) 



(where we have used 7 = 5/3 in the last form). This 
corresponds to a shock temperature of order 10 K for 
the forward shock into the interstellar medium (as op- 
posed to ~ 10 10 ~ n K for the wind shock). Since the 
outer shock and the contact discontinuity are very close 
together when energy-driven flow starts (see Fig. [1} this 
means that the outer shock is always at 



7+1 



R. 



0) 



The outflow rate of shocked interstellar gas is 
dM(R out ) _ (7 + l)/ g a 2 



dt G 
Assuming M = M a , the wind outflow rate is 

4/ c m(j 4 



R. (10) 



M w = mM Edd = 



tjcG 



(11) 



We can now define a mass-loading factor for the outflow, 
which is the ratio of the mass flow rate in the shocked 
ISM to that in the wind: 



_ A/out = V(l + 1) f g Rc 
■ M w I''' fc* 2 ' 

Then the mass outflow rate is 

7,(7 + 1) fg Rc 



(12) 



M out = f L Ah 



Mem- (13) 



4 f c a* 

If the AGN is still radiating at a luminosity close to Ed- 
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dington, we have R = v e , and using © gives 



and 



V 3ct J 



~ 460er, 



-4/3 
200 



jl/3 



3700cr^o / 1/3 M yr" 1 



(14) 
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for typical parameters, f s — f c and 7 = 5/3. If the 
central quasar is no longer active, the mass outflow rate 
evidently declines as R/v e times this expression, with R 
given by (J7J . 

It is easy to check from ([SJ [T5)) that the approximate 
equality 



1 • n 1 • 

-M w vd ~ -M, 



out "out ■ 



(16) 



holds, i.e. most of the wind kinetic energy ultimately 
goes into the mechanical energy of the outflow, as ex- 
pected f or ene r gy dri ving. One can show from the equa- 
tions in iKmel (pOOl that if the quasar is still active 
and accreting close to its Eddington rate, the shocked 
wind retains 1/3 of the total incident wind kinetic energy 
M w wd/2, giving 2/3 to the swept-up gas. The energy re- 
tained in the wind and the swept-up gas have potentially 
observable emission signatures (see Discussion). 

Equation (fTB"|) means that the swept-up gas must have 
a momentum rate greater than the Eddington value 
-^Edd/c, since we can rewrite it as 



(17) 



2M W 
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where the P's are the momentum fluxes. With P w = 
^Edd/c, we have 

A _ p Mou±\ L Edd i/ 2 -2/3,1/6 -^Edd 

out " w ) = ~ /l ~ 200 ~r 

(18) 

where /l is the mass loading factor of the outflow. The 

1/2 

factor fjj ~ 20 is the reason why observations consis- 
tently show Mout^out > L Edd /c. 

6. DISCUSSION 

We have shown that large-scale outflows driven by 
wide-angle AGN winds should have typical velocities 
v ut ~ 1000 — 1500 kms -1 and mass flow rates up to 
M out ~ 4000 M Q yr" 1 (eqns |U [TSJ) if the central quasar 
is still active, with lower values if it has become fainter. 
Our equations (0 [TBI directly relate the outflow ve- 
locities and mass rates to the properties of the host 



galaxy. The outflows should have mechanical luminosi- 
ties E out ~ (?7/2)LEdd ~ 0.05LEdd, but (scalar) momen- 
tum rates P ou t ~ 20iEdd/c- These predictions agree well 
with observations (see Tables Q] and [3]). We conclude that 
AGN outflows may well be what sweeps galaxies clear of 
gas. 

Our picture predicts several other features that may 
aid in interpreting observations. It suggests that the 
molecular outflows come from clumps of cool gas embed- 
ded in the outflowing shocked ISM. They are entrained 
by the advancing outer shock front and persist for a long 
time. We note that this shock front is Rayleigh-Taylor 
stable since interstellar gas is compressed here. Fur- 
ther, the temperature of the s hocked ISM is in the righ t 
range for thermal instability (|McKee fc Ostrikerl 19771) . 
and Richtmyer-Meshkov instabilities ()Kane et al. 1999) 
induced by the forward shock mean that new cold clumps 
may form in the outflow behind it. Simulations by Zubo- 
vas & Nayakshin (in preparation) show that up to 10% 
of the total mass in the outflow may be locked up in this 
cold phase. This agrees with the conversion factor of 
~ 10% used in the papers cited in Table[T]to estimate the 
total mass outflow rate from observed molecular species. 
Our model therefore predicts both the total mass out- 
flow rates (eqn ll5j) and the observational signatures used 
to estimate them, in good agreement with observation 
(Table p. 

The inner wind shock presumably accelerates cosmic 
ray particles, and gamma rays result when these hit the 
colder ISM and shocked wind. The outflows are therefore 
directly comparable with the gamma-ray emit ting bub- 
bles i n our Galaxy recently discovered by Fermi (|Su et al.l 
2010). One can interpret these as relics of the Milky 
Way's last quasar outburst about 6 Myr ago by noting 
that the greater density of the Galactic plane must pinch 
a quasi-spherical quasar outfl ow into a bipolar shape 
(|Su et all 120101 IZubovas et al.ll20ll . The gamma-ray 
emission from distant galaxies discussed here should be 
intrinsically stronger than in the Milky Way, but the long 
integration time required to detect the Galactic bubbles 
means that these outflows may be undetectable with cur- 
rent instruments. 

Perhaps more promisingly, these cosmic ray electrons 
cool and emit synchrotron radiation in the radio band. 
This radiation may be observable and so it would be 
interesting to check whether there are kpc or sub-kpc 
scale radio bubbles associated with the outflows. 
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